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SUMMARY 


The integral form of the velocity potential and pressure distribu- 
tion for a wing with supersonic leading edges and subsonic side edges 
in supersonic flow is derived herein for a wing undergoing any arbi- 
trary time-dependent deformations. The expressions are simplified by 
assuming harmonic deformations and then expanding the integrand of the 
velocity potential to the third power of frequency. The special case 
is treated for which the side edge is parallel to the free stream and 
the oscillations are such that the amplitude of wing distortion can be 
represented by a polynomial of any desired degree in the span coordinate 
and third degree in the chordwise coordinate. 

The equations are further reduced to the special cases of a rigid 
wing oscillating in pitch and translation and of a rigid wing in a sinus 
oidal gust, the results of which are presented in an appendix. Sample 
calculations are made for the total lift on a delta and rectangular wing 
and the results are presented in a table where a comparison is made with 
the exact values from linearized potential- flow theory. 


INTRODUCTION 


Time-dependent aerodynamic forces have been a subject of continuing 
theoretical development for many years. Most effort has been directed 
toward methods of predicting air forces due to simple harmonic motion 
since these methods can be applied directly to aircraft flutter problems 
With suitable operations these harmonically varying forces, which were 
developed for application to flutter, can be used in the harmonic anal- 
ysis of airplane response to continuous atmospheric turbulence. Thus, 
the accumulated knowledge of unsteady air forces due to harmonic motion 
of wings at various speeds may be applied to both flutter and response 
to turbulence. 
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The lift and moment for rigid wings of various planforms under- 
going harmonic oscillations have been derived. (For example, see refs. 1 
to 10.) The lift and moment on certain rigid restrained wings subjected 
to continuous sinusoidal gusts (or turbulence) have been presented in 
references 11 to 13* The results of references 1 to 13 have been com- 
piled in reference 14 together with the unsteady air forces for addi- 
tional planforms. A more complete list of references is given in the 
bibliography of reference 15. 

The aerodynamic forces for application to nonrigid or deforming 
wings are available for special cases. For example, if the distorted 
shape of the wing can be represented by a quadratic equation in the chord- 
wise and spanwise coordinates, references 1 6 to 19 may be used in the 
supersonic speed range. In reference 16, the velocity potential for a 
triangular wing with subsonic leading edges undergoing general second- 
degree forms of harmonic distortion in both the spanwise and chordwise 
coordinates is presented. The velocity potential therein is expanded 
to the third power of the oscillation frequency in order to obtain the 
forces and moments. Reference If is an extension of reference 16 wherein 
a higher degree of wing distortion is considered and the velocity poten- 
tial is expanded to the fifth power of the frequency. In reference l8, 
the generalized forces for a harmonically oscillating rectangular wing 
are given. The downwash distribution is assumed to be a general poly- 
nomial in the spanwise and chordwise coordinates. In reference 19 a 
strip theory technique is used to obtain the generalized forces on a 
delta wing with supersonic leading edges. This procedure gives the 
exact pressure distribution for arbitrary chordwise variation of dis- 
placements and, at most, linear variation of displacements in the span 
direction. 

In the present paper an integral expression is given for the pres- 
sure distribution on a wing with swept supersonic leading edges and 
arbitrarily swept subsonic side edges with an arbitrary time-dependent 
downwash distribution. The trailing edge is also arbitrary but must be 
supersonic at all points. The expression is simplified by considering 
the special case for a wing undergoing harmonic motion with side edge 
parallel to free stream. The deformed shape of the wing is represented 
by a polynomial of any desired degree in the span direction and third 
degree in the chord direction. The aerodynamic forces are obtained by 
expanding the equations to the third power of frequency. Reduction of 
the equations for application to a rigid wing oscillating in pitch and 
translation and to a rigid wing in a sinusoidal gust is presented in an 
appendix. 
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SYMBOLS 


a 

As 


b 

h 

Is 

Jn( x ) 

k 


K n 

X(x,y) 


speed of sound 

normalizing factors used in equation (46) to define displace- 
ment of wing 

wing span 

sinking velocity of wing 
quantity defined by equation (D3) 

Bessel function of first kind 

reduced frequency 

normalizing factor (see eq. (42)) 

lift distribution due to downwash, 



n=0 s=0 



quantity defined by equation (47) 




Lfc = Lk(-A) 
M 

n,s,r 

Ps 

Ap 

A p/q 


’ll 

Mach number 
integers 

quantity defined by equation (DIO) 
local pressure difference 

amplitude of pressure coefficient due to downwash, 
V = K n e^ ot x n B(y-ri) 
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q dynamic pressure, pV 2 /2 

Qs quantity defined by equation (D 4 ) 

r O-» r l> r 2 quantities defined by equations (29) 

r 3> r 4 quantities defined by equations ( 4 l) 

R s quantities defined by equation (D9) 

R(y) quantity defined by equation (37) 

t',t,t2 time 

tp transformed time (see eqs. (5)) 

u - P|y - t] | 

U = P(y + tj) 

V free-stream velocity 

w vertical velocity on surface of wing, positive up 

w 0 ,wi amplitude of vertical velocity associated with the Dirac 

delta function 

x '>y'> z ' rectangular coordinates fixed to wing 
transformed coordinates (see eqs. (5)) 

x ,y>z rectangular coordinate system fixed to apex of wing 

a angle of attack 



8(x) Dirac delta function 

1 'lS T b T l2_ position on y', y, and y^ axes where downwash is applied, 
respectively 

A slope of leading edge of wing 
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A 

P 

PM 

(£> 

- M 2 o> 
0 ) = — - 

Vp 2 

^ 1^2 


(3 times slope of side edge of wing 
density 

velocity potentials 
circular frequency 


coordinates 
dummy variable 


ANALYSIS 


Introductory Remarks 

As a first step in the analysis, an integral expression is developed 
for the velocity potential associated with a downwash strip of Dirac 
delta form on a wing in supersonic flow with swept supersonic leading and 
trailing edges and subsonic side edges. By superimposing these strips 
over the wing planform, a general expression for the generalized forces 
is derived for any arbitrary time -dependent downwash distribution. These 
equations, although complicated, can be programed on the modern-day high- 
speed digital computers. 

As a next step in the analysis, a simplification is made of the 
above-mentioned expressions by assuming simple harmonic motion. The 
velocity potential and pressure coefficients associated with the harmoni- 
cally oscillating strip are then presented. These expressions are fur- 
ther simplified so as to pertain to the special case where the side edge 
is parallel to the free stream. The pressure coefficients are then 
expanded in powers of frequency and by superimposing the downwash strips 
over the wing planform the pressure distribution for various wing dis- 
tortions is obtained. 

The method used is that of Gardner (ref. 20) which reduces the non- 
stead y finite-wing problem to two "steady" finite-wing problems. Without 
deriving the method, its essential points are given herein. 


Velocity Potential and Pressure Coefficients Associated 
With a Downwash Strip of Dirac Delta Form 


The linearized boundary- value problem .- The differential equation 
of the propagation of disturbances that must be satisfied by the velocity 
potential is (when referred to a moving coordinate system x’^y'^z 1 ) 


\bt f cht 7 


= 

dx' 2 dy' 2 dz' 2 


( 1 ) 


where 


0 s 0(x , ,y , ,z , ,t ') 

The boundary conditions that must be satisfied by the velocity poten- 
tial are 

0(x',y , ,O,t') = 0 (2) 

ahead of the wing 


(M-) =v ^-L + ^: = w(x',y-,t-) 

' ' z '-£> Sx ' St 1 


( 3 ) 


on the wing. 

The wing planform for which the velocity potentials and pressure 
distributions are to be obtained is shown in figure 1. The numbered 
regions on this figure will be discussed later in the paper. A conven- 
ient planform to analyze is shown in sketch 1 where by means of various 
transformations and superposition techniques the results can be applied 
to the planform shown in figure 1. The downwash associated with the 
planform in sketch 1 is assumed to be 

w(x',y',t') «w 0 (x',t') 8(y '-rj') (4) 

where S(y') is the Dirac delta and is defined as 
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and 


/ OO 

F(y') 5(y'-y) dy' = F(y) 

.00 


and the slope of the side edge is expressed as A/p for convenience. 



Although this dovnwash distribution may appear to be physically unreason- 
able, it will be shown in appendix A that the results utilizing it reduce 
to known functions. 


Transformation of the boundary- value problem .- With a transforma- 
tion similar to that employed in reference 20 , 


_ Ay 1 + X '/P 

1_ jr- 


y i = 


A 2 

Ax'/p + y' 


( 


1 - A c 


Mx 1 - p 2 at ' 




t l " 


z, = z’ 


ni = n’ 


( 5 ) 


j 
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equations (l) to (ij-) become 


ahead of the wing 


^ x l x l ~ " ^ z l z l ~ ^itq ° 


0=0 


= W ( X l ,y l ,t l) 


( 6 ) 

( x i < Ayi) (7) 

(y x > o) (8) 


on the wing 


r 


w ( x l^l^i) = w 0 / 


p( x i - in) _i_ 


^1 - A 2 


Pa 


M ( x i - A^i) 


i 


1 - A c 


^1 - 


111 - A 


- 


w( x l,yi,tj = v 1 (x 1> y 1 ,t 1 J 5/ 


'1 - ^1 
(/l - A 2 


7l 


(9a) 

(9b) 


where 


e( x i " ^ 1 ) 1 

"N 

Mfxi-AyJ J 

1 l/l - ^ 

1 

Tn 

> 

ro 


> (9c) 


The wing in the transformed coordinate system is shown in sketch 2 , 
where the shaded area of sketch 1 transforms into the shaded area of 
sketch 2 , and the side edge and leading edge become y-^ = 0 and 

xp = Ay^, respectively. 



c 



X 1 


Sketch 2 



and where 

x (o^x 1 ,y 1 ,t 1 ) = ( lJ 0 

and w ( x i^ 3 _^i) defined by equation (9b)* It can be seen that, by 

adding equations (10) and (ll), the resulting equation has the same form 
as equation (6), Similarly, equations (12) and (7) have the same form 
and equations ( 13 ) and (it) reduce to equation (8) as approaches 0. 

Consequently, the velocity potential is found by setting £q = 0 such 
that 
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0( x i> y i> z i’ t i) - 'i'( 0 > x 1 >y 1 > 


(15) 


Differentiating equation (ll) with respect to z^ and defining a 
new potential function X such that 


5z n 


= X(6 1 ,x 1 ,y 1 ,t 1 ) 


(16) 


z^->0 


equation (ll) becomes 


X x x - x t t - X. . =0 

x l x l t i t i Slip 


(17) 


with the boundary conditions 


0 

II 

X 

(18) 

ahead of the wing and 


x = w ( x i,yi,ti) 

(19) 

^ or £]_ = 0. This will be considered the first boundary-value 
The second boundary-value problem consists of equations (10), 
and (13) and is restated here for convenience 

problem. 

(12), 

**1*1 " ty i y i " ' kz i z i = 0 

(20) 

\|r = 0 

(21) 

ahead of wing and 


( tz i) z =0 = 

(22) 


for y > 0. It can be seen that the solution of equations (17) to (19) 
becomes a boundary condition for equations (20) to (22). 

Solution for X-function .- The boundary -value problem defined by 
equations (17), (l8), and (19) is similar to the steady two-dimensional 
supersonic "wing" problem in x-^t-^^ space, where x 1 is in the down- 
stream direction, t ± in the span direction, and f is normal to the 
wing as implied in sketch J. The downwash, as can be seen from 
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equations (9a) and ( 19 ), is concentrated along the line 



Sketch 3 

The curve x i = x 2 “ ^(t 2 - t-jJ + represents the intersection of 

the plane = 0 with the characteristic forecone emanating from the 
point 

The solution for X can now be written in terms of simple sources 



where Sq is the hatched region indicated in sketch 3* It might be 

noted that the expression for X (eq. ( 23 )) differs from the classic 
potential given for the two-dimensional steady-flow problem by the partial 
derivative d/c^. This is due to the fact that in the two-dimensional 


problem the vertical velocity, which is analogous to 




is 


specified on the wing, whereas, in this problem, the potential 
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(^)^l ~ ^ i s specified on the wing. Substituting for 

from equation (9b) into equation ( 23 ) and introducing the appropriate 
limits yields 


r 


1 a 
* 




dtc 






! y — AX 

2 - 

ifl - A 2 

l/( X l - x 2 f * (h - t 2 ) 2 - ?1 2 


dx p (24) 


By integrating equation (2^) with respect to x ^ and then making the 


substitution ^2 “ ^1 + 


. x i - 


yi - v 1 - a 

A 


- £^cos 0, keeping in 


mind that the value of the integrand is concentrated along the line 

\2 


x i = x( y i ■ ^ ~ a2 ) 

there is obtained 


between tp = t^ ± 




X = - 


VI - Isr d 




*1 


^ - a2 ^ if yi - nji - a 2 ] 2 

— j + y\ x 1 ~ j - h cos 

( A ( x i - *1) + nl - a 2 > yi > -^= 


d8 


(25a) 


and X = 0 for all other regions. When w-^ is rewritten in terms of 
w 0 ( ec l- (9c)), equation ( 25 a) becomes 



L - A c - t 1i 


■(*! " *l) - A > y l 


and X = 0 for all other regions. 

Solution for the velocity potential 0( x p>yp>°p> tp) • - The boundary - 

value problem defined by equations (20), (21), and (22) is also similar 
to steady two-dimensional supersonic wing problem in the £p>yp> z p space ; 

where §p is in the upstream direction, yp in the span direction, Zp 
is normal to the £p,yp pis- 116 as impli e( i in sketch 4. 


1|1 + A 


Sketch 4 
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The solution for \|r can now be written in terms of simple sources 
as 



and by means of equations (15) and (26) 



where Sj_ is the crosshatched region indicated in sketch 4 . The cross- 

hatched and hatched regions in sketch 4 are the regions where the poten- 
tial X(x 1 ,y 1 , 0 as dictated by the conditions imposed by equa- 

tion (25). Integration over region S]_ will yield results for pure 
supersonic flow, whereas integration over region S 2 (sketch 4 (b)) will 
yield results which contain the effects of the subsonic side and must be 

taken into account when 0 < y^ < x-^ - 



By substituting the appropriate limits into equation ( 2 j), 
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where 


r 0 = 


A ( x i + y i) + \r ” aC 

1 + A 




r l = 


r 2 “ 


niV 1 - A - y 2 


X 1 


A ( x i ~ y i) + ^iv 1 ■ AC 

_ i + A 


(29) 


If the potentials in equations (28) are transformed back to the true 
x'^y'jt' coordinates, the regions where these potentials exist are 
given below and are also shown in sketch 5 


For tj* > 0: 

0 = 0 1 (x',y',t ') 




tj ' - — < y ' < t) ' + — , 0<x’< — 

\ P P 1 - A 


0=stU',y,t') |^-v(H 4 ) 


< y ' < q ' + — , x ’ > ^ 


*1* 

- A_ 


> ( 50 a) 


U -J J 


and 0=0 for all other regions. 
For r\ < 0: 


= 0 2 (x',y',t') l- — <y’ <-J A 


Ax' ^ y x' 1 + A , , Pri T \ 

+ — ■ — n\ x f > - i-l- 


A 


( 50 b) 


and 0=0 for all other regions. 


[g 


0=0 


Wing side edge 

y =-^ 


0 = 0c 


'<///////// 

<^W(x',y',t'X 


(a) T] ' > 0 



y' = — - n' 

X ' p 




/ / / / / / / 


Wing side edge 
. Ax r x 1 

y 


(b) T)' < 0 


Sketch 5 

The pressure difference in terms of the physical x'^y'^t' coordinates 
is given by 


Ap ■ 2p W + v & 


Subsequently, a coordinate system is employed where the origin is 
situated on the apex of a sweptback wing. The x-coordinate is parallel 
to the free-stream direction and the y-coordinate, in the span direc- 
tion. The coordinate system will also have the additional character- 
istic that the Dirac delta strip will have its origin on the leading 
edge of the wing rather than on the y'-axis as indicated in sketch 5(a) 
If it is assumed that the potential functions given by equations (28) 
and the pressure coefficient as given by equation (31) have been trans- 
formed to this new x,y,t coordinate system, the lift distribution due 
to arbitrary downwash distribution can be written as 

Hj , pk 

A Pl (x,y,t;n) Ap 2 (x,y,t; n ) ^ 


7 ( v 
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where Ap^q and Ap^q are the pressure coefficients associated with 
the velocity potential 0^ and 0 2 , respectively, and f(q) is a non- 
dimensional spanwise downwash weighting function defined as 

v(x,y,t) = w Q (x,t) f(y) 

A method for obtaining the limits of t) will be discussed in a subse- 
quent section; however, an indication of their significance is shown in 
sketch 6 where r^, T) ., and are particular values of t) s 

y>t\ 


x 

Sketch 6 

It might be noted that, in order to obtain Ap/q for use in equation ( 32 ), 
a triple integration is involved: one as indicated in equation (25b) and 

two more as indicated in equations (28). Therefore, in order to obtain 
the loading due to an arbitrary downwash distribution, four integrations 
are required. To obtain the generalized forces another integration is 
required, and if a transient phenomenon is present, a time-superposition 
integral is required. However, with the modern high-speed computers it 
does not seem unreasonable to undertake the integration of a quintuple 
or even a sextuple integral. 

The remainder of this paper will deal with the evaluation of equa- 
tions (25b) and (28) in order to obtain the lift distribution for a wing 
oscillating in simple harmonic motion with a polynomial downwash distri- 
bution in the chordwise and spanwise direction. The potential will first 
be derived for the oscillating strip, whose downwash can be represented 

by 
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w ( x Sy % t ' ) = v 0 (x',t') 6(y ' - r| ' ) 

= w 0 (x ' ) e^' &(y ' - T, 

The potential is then derived for this strip and a frequency expansion 
performed to obtain the corresponding pressure distribution. These 
expressions are further simplified for the special case where the side 
edge is parallel to the free stream and the downwash is assumed to vary 

as x ,n in the free-stream direction. Superposition techniques are 
then used to obtain the loading distribution for various wing distortions. 



Velocity Potential and Pressure Coefficients 


for Simple Harmonic Motion 


Solution for the potentials X and 0 .- By assuming simple har- 
monic motion, the downwash distribution as given by equations (33) becomes 


w ( x l^l^l) 




io) 

M ( x i- Ay i) . 


p( x l - Ay 1 ) 

pa 

e 

1 — 1 

_ yi-A 2 

0 

^1 - A 2 






( 3*0 


when transformed to the x-^y-^t^ coordinate system by means of equa- 
tions ( 5 )* Comparing equations (3^) an( l ( 9 ) gives 


w o( x i^i) = w o( x l) 



Substituting equation (35) into equation (25b) gives 


( 35 ) 



(36a) 
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and X = 0 for all other regions or 


X , -R( yi )— Jo 


CD 

3a I 


» X 1 - 


y i - V 


1 - A c 


A 


+ n\ll - A 2 > Y, > ^ 


a( x i - Ijl) + VjV 1 - A > y± 


l 


1 - A 


( 56 b) 


and X = 0 for all other regions where 



If the indicated differentiation is performed and it is kept in mind that 

Ax x + T]j\jl - A 2 - y-L 

there is a discontinuity in X along the line = T > 

equation ( 36 b) becomes 



The evaluation of potentials 0-j. an<i 0 2 ( e< l s - ( 2 ^)) together with equa- 
tion ( 58 ) is presented in appendix A and the final results for the veloc- 
ity potential (eqs. (A29) and (A30)) in the true x'^y'^t' coordinate 
system are as follows: 


0 1 ( x ',y , , t ') = - 


1 icDt T 

— e 

it 



v 0 (x'-3£) 


-iccMl/p 


a 

cos T 





(39) 
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(^ 0 ) 


where the regions where these potentials exist are given by equations (30a) 
and (30b) and are shown in sketch 5 and 



Frequency expansion of velocity potential and pressure coefficients . - 
An analytic evaluation of equations (39) nnd (4o) does not seem possible 
at the present; therefore, a downwash distribution 


w ( x '>y'>t') = K n x ,n e Jxot &(y'-T] ' ) (42) 

is chosen, and a frequency expansion of the integrands of equations (39) 
and (40) is performed before integration. It might be noted that com- 
parison of equations (42) and (33) indicates that 


v 0 (x') = K n x' n 


(^ 3 ) 


21 


where K n is a normalizing factor having the dimensions of velocity- 
divided by the nth power of a length. 


The results of the above-mentioned expansion are presented in 
appendix B. Equations (B5) to (B22) represent the pressure distribu- 
tion for a wing with a supersonic leading edge and subsonic side edge 


r) 


in supersonic flow. The downwash distribution on this wing 


('•-T 

is given by equation (4-2). By superposition techniques, the pressure 
distribution over the entire wing for any harmonic deformations can be 
obtained; however, the amount of work and time involved becomes very 
lengthy. Therefore, the pressure coefficients are derived only for the 
special case where the side edge was parallel to the free stream. These 
coefficients are presented in appendix C as equations (Cll) to (C2l) 
for values of n = 0, 1, 2, and 3- The values of n = 0, 1, 2, and 3 
represent a chordwise strip dy' of the wing at y' = tj 1 undergoing 
translation, pitching, parabolic bending, and cubic bending, respectively. 


In most analyses it is desirable to have the origin on the center 
line of the wing. Therefore, in the section to follow a coordinate sys- 
tem is chosen so that the origin is at the apex of a sweptback wing. By 
assigning a given spanwise variation of deflection and integrating over 
the appropriate region of the wing, the pressure distribution can be 
obtained for any spanwise variation of deformations and up to a cubic in 
chordwise variation of deformations. 


Loading Coefficients for Polynomial Downwash Distribution 

Transformation to a coordinate system fixed on apex of wing . - A 
sketch of the wing together with the new coordinate system, fixed on the 
apex of the wing, is shown in sketch 7 where A is the slope of the 
leading edge and b is the wing span. Two sketches are needed depending 
on whether rj * is greater than or less than b/2. 



y 


y' 



x' 


(a) q’ < b/2 


Sketch 7 


(b) q 1 > b/2 
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Inspection of sketch 7 shows the appropriate transformation to be 



The resulting figure showing the wing, new coordinate system, and the 
regions where the appropriate pressure coefficients apply is shown in 
sketch 8: 



Sketch 8 

and the downwash equation becomes 

w(x,y,t) = Kne^^x - hij 5(y-ri) (45) 


General loading coefficients .- If the wing is now divided into 
regions according to characteristic Mach wave reflections (see sketch 8) 
and the downwash is given by 



w(x,y,t) = K n e 
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then the pressure at any point will consist of terms of the form 




Z Ap . 

Vl S — ( x -^ y + p n +|J dr) (47) 

hi 

where tj^ and tj are to be chosen according to the region in which 

the pressure coefficient is to be calculated and Ap^/q values are given 

in appendix C. As an example, consider the pressure at point x,y indicated 
in sketch 9 



Sketch 9 

from which it can be seen for this particular case 


% 


= A 


Py - * 

Ap + 1 


|j( b - y) - 3 

* A — 5£Ti — 


The pressure at point x,y can now he written as 
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i(x, y) = 


I 






K - 


y 


+ 


b 

2' 





d n 



(48) 


The last term has the form of equation (47) and can be replaced by 
The first term 


b/2 

o 


I 


As 


PA 


~P(b-y)-x 


PA+1 


(3A+1 




when A is replaced by (-A). 


This expression is now in the form of 


equation (47) and can be represented by 


nj (-?0 

where the bar (-) 

ti ± (-A) 


is used to indicate that A roust be changed to (—A) before inserting 
the limits. This definition was adopted so as to utilize the definition 
of given by equation (47) and thus eliminate the derivation of a 


new set of equations. Performing the same operation on the second term 
of equation (48) results in the following expression for the loading 
coefficients at x,y. 
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I 

1 

0 

l(x,y) = L x 

+ L 2 

tI-l(-A) 

+ L 2 


b/2 


( 49 ) 


The Lj_ and functions are derived in appendix D for 

n = 0, 1, 2, and 3* The limits are left arbitrary so that the pressure 
coefficients at various positions on the wing (other than the position 
considered above) could be determined. 

Loading coefficients for separate regions on the wing .- If the wing 
is now divided into the regions indicated in figure 1, the loading coef- 
ficients for each region ^ Zjp etc.) can be shown to be 


l T - L~i 


% 


(50) 


Z II ~ L 1 


+ L„ 


^0 


b/2 


(51) 



^2 

^0 

b/2 

^0 

Z III - L 1 

- h 

+ L 2 

= 2 II + L 1 

- 1 ! 


0 

0 


0 


^0 


0 


(52) 


2 IV “ L 1 


r\ 2 

0 

+ L 2 

+ L 2 

^0 

^2 


b/2 


(53) 


h ~ L i 


o 


+ L, 


*0 


“ 2 I " 2 II + 2 m 


(54) 


0 
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0 n 2 ^3 b / 2 

Zyj = 

^3 0 -b/2 r) 2 

^2 n 5 o b/2 

l VII = L 1 + L 2 + L 2 + L 2 = l Yl + l iv 

n 5 -b/2 Tjg 0 

where 


(57) 


and the values of and Lg are given in appendix D as equations (Dll) 

to (Dl8) and, as pointed out previously, the bar (-) indicates that A 
must be changed to -A before substituting in the limits. 

RESULTS AND DISCUSSION 



(55) 


- l IT1 (56) 


The integral expressions for the velocity potential and pressure 
coefficients associated with a wing with swept supersonic leading edges 
and arbitrarily swept subsonic side edges, deforming in any general time- 
dependent manner, are derived herein. The expressions are simplified by 
first assuming harmonic motion and then expanding to the third power of 
frequency. The special case is then treated for which (l) the side edge 
is parallel to free stream and (2) the oscillations are such that the 
distortion of the wing can be represented by a polynomial of any desired 
degree in the span coordinate and third degree in the chordwise coordinate. 
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Calculations are made to obtain the total lift coefficients for two 
wings - a 50 ° delta wing and a rectangular wing of aspect ratio 0.8 - both 
flying at a Mach number of J.O. The wings are assumed to be subjected 
to continuous sinusoidal gusts and to harmonic sinking oscillations. No 
spanwise variation in downwash is considered. This analysis is presented 
in appendix E and the final results are tabulated in table I . As can be 
seen. in table I, the results are in good agreement with those obtained 
by using reference 14 for reduced frequencies at least as high as those 
indicated in the table. It might be noted that the results for the wing 
in a sinusoidal gust are not as good as those for a wing undergoing har- 
monic sinking oscillations. It is believed that this difference is due 
to the fact that the sinusoidal gust wave is approximated by a cubic in 
the chord direction, whereas the downwash for harmonic sinking oscilla- 
tions is exact. 


CONCLUDING REMARKS 


The integral expressions for the velocity potential and pressure 
coefficients associated with a wing with swept supersonic leading edges 
and arbitrarily swept subsonic side edges, deforming in any general time- 
dependent manner, are derived herein. The expressions are very compli- 
cated; however, with the modern high-speed computers it does not seem 
unreasonable to undertake such a task. As a possible check the equations 
are simplified by assuming simple harmonic motion and expanding to the 
third power of frequency. The special case is treated for which the 
side edge is parallel to the free stream and the oscillations are such 
that the distortions of the wing can be represented by a polynomial of 
any desired degree in the span coordinate and third degree in the chord- 
wise coordinate. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 1, 1962 . 
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APPENDIX A 

REDUCTION OF THE VELOCITY POTENTIAL EQUATION 
WITH APPLICATION TO A RECTANGULAR WING 


Reduction of Velocity Potential Equations 


The velocity potential as defined in the text by equations ( 28 ), 
( 37 )> and (38) will be treated in six parts as indicated below. From 
equation (38); let 




(Al) 


where 


X i( 6 i' x i>yi'ti) = r (yi) sU - 


Ax x - + qjl - A 2 


A 


^(h’Xl’yi’h) = R ( y l)^~ J 0 


(JO 

Pa 


lx l - 


y i " ^iV 1 ■ 


- I- 


> (A 2) 


'J 


By substituting equations (Al) and (A 2 ) into equations (28) and defining 
the quantities 


e i(^2^i^i^i; y 2 ) = j==^= 

p 2 2 - (n - y 2) c 

p 2 2 - (y x - y 2 y 


(A3) 


the velocity potential of equations (28) can be written as 

0l(xi,yi, o ,ti) = 0 n - jZS 12 


(Ah) 


^ lAr^O 
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02( X l ,y l ,O,t l) = ^21 + 022 
0 21 = $211 “ ^212 
022 = 0221 “ 0222 


where 




0221 


1 

It 



r 


2 



dy 2 



at 2 


(A5a) 

(A5b) 

(A5c) 

(A6) 

(AT) 

(AS) 

(A9) 

(A10) 
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(All) 


and Tq, r^, and r ^ are defined by equations ( 29 ). 

Since the area of integration does not include the li ne 

+ n-fi - a2 - y 2 

= — (which is the argument of the Dirac delta). 


0221 


= 0 


(A12) 


After the integration with respect to equations (A 6 ) and (A 8 ) may 

be represented by 




1 


pr. 


R ( 

£ 2 ] 

1 <ty 2 

jf* 1 - y 2 + ^ 

E 


2 

| " ( y l “ y 2) 2 


(A13) 


By transforming equation (A13) to the true x',y' coordinates by means 
of equations ( 5 ) and making the substitution 


y 2 



-AS + q’ + 


Ax 1 \ 

P / 


equation (AI 3 ) becomes 



(AlU) 


fl 2 - (y* - q ') 2 


(A15) 
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With the appropriate substitutions for r^ and r^ it can be shown that 
equations (A6) and (A8) become, respectively. 




In order to reduce equations (A7), (A9)> an d (All) the expression 



will be used to represent these three equations. By means of the 
substitution 






( y i - y 2) 2 


(A19) 


equation (Al8) takes the form 


32 


R (yg ) fog 




*2 


- nj 1 - A 2 ) - (y x - y 2 ) 


(I - r 2 ) 


A 7 




( y i- y 2) Jo 


a 7 


(A20) 


which, after transforming to the true x^y^t 1 coordinates by means 
of equations ( 5 ) and making the substitution given by equation (Al4), 
reduces to 


K 


1 ia>t' 
— e 
it 


_ 

w o( x '-pO e Pa 
- (y' - n ’) 2 


d£ 



03 p l03t 




(A21) 


An expression that will be needed later is that for r^ = 1 the 
second integral, in equation (A2l) becomes 



(A22) 


where 

T = fife 2 " ^ y ’ " (A25) 


By letting 7 = sin 9, the integral (A22) becomes 



sin 9 ) d 0 


(A24a) 
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which is evaluated by Watson (ref. 21, p. 37*0 and is given as 

i(l - cos t) (A24b) 

By substituting the appropriate values for r^, r^, and r^ into 

equation (A2l), it can be shown that 

r x '/P 


t = _ i ' 

led It 


-(y'-V) 


•W 0 (x'-P0 e (l - cos t) d£ 

k 2 - (y' - n ') 2 


= _ 1 loot ' 


212 « 


Del 

Wq(x'-PO e (l - cos t) 

\k 2 - (y' - V ) 2 


-(y'-n') 



Consequently, if equations (A12), (Al6), (A17), (A25), (A26), and (A27) 
are substituted into equations (A4) and (A5a), the resulting form for 
the potential is 


3 ^ 


0i ( x Sy St ') = - i e icot ’ 


■73 


ly’-n'l 


_ iopM ^ 

w 0 (x'-pO e Pa ‘ 


COS T 


{P - (y* - p ') 2 


(A29) 


0 2 ( x ',ySt') = -ie 

Jt 


5 IcuMf . 

" Pa 


1 „io)t 1 


Wq(x’-P^) e p cos t 


y'-n'j |/5 2 - (y 1 - V) 2 


d5 


-x’/p 


+ ^ e io)t * 

*0a 


kuM„ 


w q( x ' - P£ ) e Pa 


' r 4 


1(4 - s) 


0 £7 


d7 (A50) 


Application to a Rectangular Wing 

It might be noted that, as A approaches 0, equations (A29) and 
(A30) become 


^i)a=o = 


1 e icDt ' 
n 


rx 


7p 


W 0 (x'-P£) e p cos t 


iiuM 

Pa 


t 


J 


|y'-V| F 1 <y ' - ’ ,)2 


(A31) 
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(feu 


1. e io)t 1 
it 



_ lmM ^ 

w 0 (x'-pQ e pa cos T ^ 

\t, 2 - (y' - n') 2 


+ e lu)t ' 
Ttpa 




(A32) 


It is now desired to find the potential for a rectangular wing for which 

w(x',y',t') = Wq(x') e^'. Since equations (AJl) and (A32) apply only 

for a downwash strip in the free- stream direction at y' = t]', these 
potentials must be integrated with respect to ’ over the proper limits. 
From an examination of sketch 10 



Sketch 10 


the potential for the rectangular wing can be written as 



0x'<py' “ 


(A35) 


36 



(A34) 


An evaluation of the integrals in equations (A33) und (Aj4) vill yield, 
the results given by equation (15) of reference 10. 
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APPENDIX B 

FREQUENCY EXPANSION OF THE VELOCITY POTENTIAL 
AND PRESSURE COEFFICIENTS 


Substitution of the downwash wq(x') = K n x ,n (eq. (L 3 )) into the 
expression for the velocity potential 0^ (eq- (39)) and expanding the 

integrand in powers of frequency, and integrating yields the following 
results (where the primes have now been omitted from x, y, t, and q) 


0 1 U,y,t) = - | e“ ofc K n ^(-icu) ] 


hpNp 


(Bl) 


p=0 


vhere 


h P = 




l?/S n s+ r 

Y 1 V (-D S nl(p|y-l|) P y (-l) m r 

r 4o M 2r (2r)l(p - 2r)l A, sl(n-s)! Zj ( r . m ) 

-1 X 


(r - m) !ml 


> (B2) 


\iosh 


N p = x n - s 


P I y-Ti i 


s+p-2m, , 

cosh 0 d 0 


0 


J 


— M^cjd 

OJ = — 

Vp^ 


(B3) 


and [p/2] is the integer part of p/2. The pressure coefficient asso- 
ciated with is 


ffi = ^ k M + v 

9 q y2\c)t 5x / 


(B4) 


58 


or 


A Pl 4 


S/ K " - Up' 

p=0 L M j 


and 




M-Sl X 


s+p-2m 


N ' = 5? = !Lj4ji » + 

P dx x P 


p-2m 


- M, 


x 2 - 3 2 (y - T ]Y 


(n £ 0) 


N ' = 
P 


Ply - nl 


l/* 2 - P 2 (y - I) 2 

Similarly, equation (40) can be expanded to give 

-^K n 


(n = 0) 

J 


Ap 


p=0 


(B5) 


> (b6) 


(B7) 


^.a^K.^-a-y 


S m-q F q + g m-q F q -itu )g m -q F q 


m— 0 q=0 


where the prime indicates differentiation with respect to x, and 


AP 2 ~ AP21 AP22 


(B8) 


(B9) 


>cosh 


-1 


N p = x n-s 


ha.] 


(l+A)(y-q) 


s+p-2m 

cosh 0 d0 


(BIO) 


0 


39 


2p n +2| 

(-♦£) 

m-q 

1 i 

K* + f ) 

4M^( -A) n+ ‘ L (in - 

■ q) ! I 

h- ' 

1 

> 

ro 


(fill) 


2 n 

Fq = ^ b q-2k ^k ^ c s G q-2k,s,k ( B12 ) 

k=0 s=0 


G <l,s,k J 

r ’i q+s + i 

H 2 f k (l-0 <i|i 

r 5 

(B13) 

1 - A 

r S = 

5 1 + A 

q(l + A) - A^y - |j 

(Blk) 

t 

II 

“S 

(B15) 

c s = 

/ ,is . n-s 
(-1) n! T) 

s ! (n - s ) I 

(B16) 

2k 

■ ID 

(BIT) 


fo(^) = 1 

(B18) 


-§(- 8 -?) 

(B19) 

^ - 2131 

/ 2 2' 
-lz k - 2 7UZ 2 * 7 M 

V 5 5 > 

) (B20) 
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(B21) 


(B22) 


With the limited number of f k (p) functions given (eqs. (Bl8) to (B20)) 
equation (b 8) can only be expanded to the fourth power of frequency. 



APPENDIX C 


FREQUENCY EXPANSION OF PRESSURE COEFFICIENTS 


FOR SIDE EDGE PARALLEL TO FREE STREAM 


Since equations (B5) and (B7) of appendix B are independent of the 
slope A attention need only be given to equation (B8). However, since 

Ap 22 

singularities exist in equation (B8), was evaluated by using 

1 A=0 

the second term in equation (AJ2) of appendix A. The resulting equations 
are (where the primes have been omitted from x, y, t, and tj) 


f 22 

I A=0 




,q,k s,q,k 


where the prime indicates differentiation with respect to x and 


CL 


L L 

k=0 s=0 


c s d k K s ,q-2k,k 


K s,q,k x 


3(y +, n) 


ct^ +S -3^( ct) da 


^ p! 


(-D nl 
(n - s ) 1 s! 


-m 

\2M 


■a 0 (tr) = 1 





(c8) 


^ (o> “ 5T31 fi 4 - ^ z i £ + i ^) < c ?> 

z x 2 = o 2 - p 2 (y - T|) 2 (CIO) 

p 2 L J 

Evaluating the above expressions for values of n = 0, 1, 2 } 3 and 
retaining only terms to the third power of frequency, the following 
results are obtained. In the following equations u = p|y - rj| and 
u = p(y + Tj). 

For n = 0: 




^3 






q. 


kK ± 

W 


+ 


C 1£ t 
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- 12M 2 + j)xu 2 + i(5M 4 - 6M 2 + l)uu 2 > 
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^ { % 0^(« - S) 2 + P^(x - a) 2 [x ♦ ( 5 M 2 - l)u 


For n - 3- 

Ap-L 

q 


4K, 


^x 2 + ^cosh’ 1 ^ | x'Jx 2 - » 2 ] - ^ ^|2p 2 x 5 + (?M 2 - l^xu^ cosh' 1 ~ 


+ ijx 2 - u 2 (lTM 2 - ll)x 2 + l^M 2 - l)u 2 | > + -^-5 <-u 2 p63 2 x 2 + }(5M 2 - 3)u 2 


i6?r 


cosh 


•1 x 


^x 2 - u 2 2 (m 2 -3 )x 2 + ( 149 M 2 - 39)u 2 j\ - Jl53 2 xu 2 J^3 2 x 2 + ^M 2 - iju^jcosh" 1 ^ 


a P2 I 


4k : 

rtV 


+ U 2 - U 2 [2 (m 4 - 6M 2 - 3)x^ - (l59M^ - 274m 2 4- 83)x 2 u 2 - ^(sm 4 - 8m 2 + l )\ i * 


— (l3x 2 + ^-Jcosh' 


IT 1 H + 20 ^7| (- 6 x + | u j - ^ <^x|p 2 x 2 + j(jM 2 - l)^cosh -1 i± + 2 P^ 

. |(jm 2 - i)xu + (^n 2 ^- jj 52 + ?( 1 > m 2 ~ 1 ) u 2 j> + ^ i |(^ 6p2x2 + 5 ^ 5m2 ” ^) u ^| cosh l 4 

+ 23^[p 2 x3 - x2a ♦ 2M 2 xu 2 + (UM 2 - 5 )xu 2 - 5? - uu 2 ]| 


i(2M 2 


iSj5 jp 2 xu 2 


1 16 


|^3 2 x 2 + ^(^M 2 - l)u 2 cosh - ^ ^ + 2p\jyr[j^ 


■fa-l 3- i+M 4 + 3M 2 - 3 2-2 


+ ^M 4 + 6m 2 - 3 ^3 _ 2 M 4 4 - 3M 2 - l a 4 _ 8M 4 - I2M 2 + 3 ^2 + 3^ - 6 m 2 + l) ^2 
8 10 6 1 6 


8M 4 - 8M 2 + 1 -2 2 8M 4 - 8M 2 + 


1 


30 


13 


Ap 22 _ ^3 

q ~ ~ ttv 


I g p^x - u) ? + pvr(x . 4 + (i» M 2 . i)ujj> 


(C19) 


(C20) 
- l)x 2 


(C21) 


(C22) 
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Equations (Cll) to (C22) now represent the pressure distribution 
on a wing (in supersonic flow) with a supersonic leading edge and a 
side edge parallel to free stream for which the downwash distribution, 
as mentioned previously (with primes omitted) is 

w(x,y,t) = K n e iayfc x n b(y-T]) (n = 0, 1, 2, and 3) 

The values of n = 0, 1, 2, and 3 represent a chordwise strip of the 
wing at y = rj undergoing translation, pitching, parabolic bending, 
and cubic bending, respectively. The wing and the regions where 2p/q 
apply are presented in sketch 11. 



x 


Sketch 11 



APPENDIX D 


EVALUATION OF EQUATION (47) 


Equation (47) as given in the text is restated here for convenience. 
The two forms are 



Examination of Ap^/q in equations (Cll), (Cl4), (C17)> and. (C20) 
indicate the presence of two particular integrals I s and Q s which 
are associated with L^_ and will be defined as 


Is = 




% 


n s d n 

x " a) " p2 ^ y " ^ 


pT1 j 


^S “ 


(y - 


x ~ T 

q ) S cosh“"*' — dT) 
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(D4) 
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where 
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— - \|(Ax - T ]) 2 - 0 2 A 2 (y - 
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'J 
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+ (2s - l)(p 2 A 2 y - Ax)l s _ 1 - A 2 (s - l)(p 2 A 2 - x 2 )l 
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= - (y - n) 


s + 1 


T ) 

3+1 1 x - 7 

cosh _± A 


P|y - n| 


(D7) 
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y 

+ i) L 


sly s - r (-D r 


n=ni 


A(s + 1 ) i—> rl(s - r) 1 r 
r=0 


(d8) 


Examination of the expressions for 6p 21 /q. and ^p 22 /q given in appen- 
dix C indicate the presence of two particular integrals R s and P s 
associated with Lp and will be defined as 


R s = 


¥ 


- n J !£ - ndn = 




s+3/2 


S + 2 


(D9) 
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?S 
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(DIO) 


By substituting equations (Cll) to (C2l) into the appropriate 
expression for Lq or Lg and using the quantities defined by equa- 
tions (Dj), (d4), (D 9 ), and (DIO), the following equations for Lq and 
Lp are obtained. 
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(D17) 
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APPENDIX E 

EXPANSION OF DOWNWASH FUNCTION 


Rigid Wing in Continuous Sinusoidal Gust 
In order to utilize equations (Dll) to (Dl8) for a wing in a sinus- 

Wt-£) 

oidal gust, it is assumed that the dovnwash w = v^e ' v ' when expanded 

to the third power of x, will adequately define the gust. Although, 
in general, this is not true, a very good approximation to the gust func- 
tion can be made if the reduced frequency k is restricted to permit 
only a third of a wavelength of the gust to be on any chordwise strip 
at any instant of time. This restriction on k is not unduly severe 
since for most analyses the useful frequency range is well within the 
limits stipulated. 


On the basis of these assumptions, the downwash can be rewritten 


as 


w = w Q e 


K^v). 


= w Q e IiLlt e 1?ocp « w Q e i4u1 '(l - iApx - P x + i fhAl) (El) 

where Ap = — By means of the following identities, 
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■Hr -K-9 -Sf-9-g 

equation (El) can be rewritten as 
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Examination of equation (E2) and equations (Dll) to (Dl8) shows 
that for a wing in a sinusoidal gust the values of A s can he defined 
as follows: 

for equations (Dll) and (D12), 


Aq = 1 

A x = -ip 



for equations (D13) and (Dl^), 


Aq = -iAp 


A-, 



_ ip-^A 


for equations (DI5) and (Dl6), 


Aq = 


2-2 

-p A 


1 2 


and for equations (D17) and (Dl8), 

_ ipV 

*<> - ~r~ 

Here the unit of length associated with was chosen as one so that 

Kq _ wq 

~ *” - 1* On the basis of these values of the problem was 

set up on the IBM 65O data processing machine to retain only the third 
power of frequency. For example, in equations (Dll) and (D12) all the 
terms in the bracket were retained for s = 0, whereas only the first 
three terms for s = 1, the first two terms for s - 2, and the first 
term for s = 3 were retained. 
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Oscillating Rigid Wing 


This abbreviated program is not restricted to repeated sinusoidal 
gusts. For example, for a, wing undergoing sinking and pitching oscilla- 
tions, the downwash can be written as 


w = e ia * 


Va 


+ h + ia^x - XQ^a 


= Ve 


iayt 


h 

ry -f _ _ 

iaiXQa 

+ v + kca/v 

Z\ 

* V 

V 

VA 7 + ~V 

a)_ 


(E3) 


where, for equations (Dll) and (D12), 

Aq = OL + | - 


h liOCUC 0 


A , - 

1 VA 


and for equations (D13) and (Dl4), 


Aq = 


ioxL 


All other values of A s are zero. Again the unit of length associated 
with K n was chosen as 1; thus, K n /v = 1. 


Application to a Delta and Rectangular Wing 

Limited forms of equations (5^); (5^); and (5*0 have been programed 
on the IBM 65 O data processing machine. The amount of information pro- 
gramed was dictated by the form of the downwash function for a wing in 
a continuous sinusoidal gust field as indicated by equation (E2). In 
order to check out the program, calculations were made for two wings, 
a 5 O 0 delta wing and an almost rectangular wing of aspect ratio 0.8, 
both flying at a Mach number of 3.0. The tangent of the leading-edge 
sweep of the almost rectangular wing was 10“6 instead of zero, since 
for A = 0 singularities arise for which no provisions were made in 
the program. The position of the points for which the pressure coeffi- 
cients were calculated was determined by means of a Gaussian distribu- 
tion formula. 

Four semispan stations were used and each region along a chordwise 
strip was divided into three stations. This procedure allowed for 
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calculating the pressure coefficients at l 8 points on the half-span 
delta wing and 21 points on the rectangular wing as indicated in 
sketches 12 and 13 . 
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Sketch 12 


Sketch 13 


This procedure for positioning the various points permitted the use of 
a fifth-degree (three point) Gaussian integrating formula for each chord- 
wise region and a seventh -degree (four point) formula spanwise to obtain 
the total lift on each wing. In table I the results are presented for 
a 50° delta wing flying at a Mach number of 3-0 in a sinusoidal gust 
field. As can be seen, the total lift coefficients are in good agree- 
ment with those obtained by using equation ( 67 ) of reference 14 for 
values of reduced frequencies at least as high as those indicated in 
the table. 

As a further check, the total lift coefficients for a 50° delta 
and an almost rectangular wing of aspect ratio 0.8, both at a Mach num- 
ber of 3 and undergoing harmonic sinking motion (a = 0 in eq. (E3)), 
are also presented in table I. The results obtained for these cases 
were in very good agreement with those obtained by using equations (66) 
and (71) of reference 14 in the same frequency range. 
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.2222 1.4093 - .0274i 1.4093 - -02551 .2667 1-0988 + .02401 1.1001 + .02331 

•3333 1.4025 - .0396i 1.4035 - . 0370i .3556 1.0986 + .03141 I.0992 + .0322i 











